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Abstract 

We calculate the imaginary parts of the isoscalar scalar and isovector electromagnetic 

form factors of the nucleon up to two-loop order in chiral perturbation theory. Particular 

attention is paid on the correct behavior of ImcTjv(i) an d ImGg M (i) at the two-pion 

threshold to = 4ml in connection with the non-relativistic 1/M-expansion. We recover 

CN ; the well-known strong enhancement near threshold originating from the nearby anoma- 

■ lous sing ularity at t c = 4m 2 ~ m*/M 2 = 3.98ml. In the case of the scalar spectral 

function Imcrjv(t) one finds a significant improvement in comparison to the lowest order 

one- loop result. Higher order 7T7r-rescattering effects are however still necessary to close a 

remaining 20%-gap to the empirical scalar spectral function. The isovector electric and 

magnetic spectral functions ImGg M (t) get additionally enhanced near threshold by 

the two-pion-loop contributions. After supplementing their two-loop results by a phe- 

' nomenological p-meson exchange term one can reproduce the empirical isovector electric 

C ■ and magnetic spectral functions fairly well. 
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^ ; PACS: 12.38.Bx, 12.39.Fe, 13.40.Gp. 
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The structure of the nucleon as probed in elastic electron-nucleon scattering is encoded 
in four form factors G S ^ V M {t) (with t the squared momentum transfer). The understanding 
of these form factors is of importance for any theory or model of the strong interactions. 
Abundant data on the electromagnetic form factors of the nucleon exist over a large range of 



momentum transfers, < \f—i < 6 GeV. Dispersion theory is a tool to interpret these data 
in a largely model- independent fashion. Each electromagnetic form factor G s ^ v M {t) can be 
written in terms of an unsubtracted dispersion relation and their imaginary parts are often 
modeled by a few vector meson poles. This procedure is based on the successful vector me- 
son dominance hypothesis which states that the (virtual) photon couples to hadrons only via 
intermediate vector mesons. However, as pointed out already in 1959 by Fulco and Frazer [1] 
such an approach is not in conformity with general constraints from unitarity and analyticity. 
In particular, the singularity structure of the triangle diagram (left graph in Fig. 1) is not re- 
spected at all when using only vector meson poles. As a consequence of a nearby logarithmic 
singularity at t c = 4ml ~ m\jM 2 (on the second Riemann- sheet) the low-mass two-pion con- 
tinuum has a very pronounced effect on the isovector electric and magnetic spectral functions 
lmG v EM {t) on the left wing of the p(770)-resonance. This becomes particularly visible in the 
isovector mean square radii of the nucleon. The effect was quantified by Hohler and Pietarinen 
[2] in their work on the nucleon electromagnetic form factors based on tttt — > NN partial wave 
amplitudes. The dispersion relation analysis of the nucleon electromagnetic form factors was 
recently refined by the Bonn-Mainz group [3, 4] accounting for new electron scattering data 
and the high-energy constraints from perturbative QCD. 

The first one-loop calculation of the isovector electromagnetic form factors has been per- 
formed in ref. [5] using the fully relativistic version of baryon chiral perturbation theory. The 
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correct analytical structure naturally emerged in this calculation and thus the strong enhance- 
ment near the two-pion threshold. In the framework of heavy baryon chiral perturbation theory 
next-to-leading order one-loop results for the isovector electric and magnetic spectral functions 
ImGg |M (l) have been given in ref.[6]. To that order a formally incorrect threshold behavior, 
not following the required p-wave phase space \m.G\ M (t) ~ (t — 4m^) 3//2 , has been obtained. 
This deviation is a consequence of the coalescence of the normal threshold t = 4m 2 and the 
anomalous threshold t c = 4m 2 — m\jM 2 in the heavy nucleon limit M — > oo. In that paper 
[6] also the isoscalar electric and magnetic spectral functions ImG^ M (t) (starting at two-loop 
order) have been computed. It has been shown that no substantial contribution arises from 
the three-pion continuum below the narrow a;(782)-resonance. 

The purpose of this paper is to present results for the isovector electromagnetic spectral 
functions up to two-loop order in chiral perturbation theory paying special attention on the 
correct threshold behavior. We will include in our discussion also the isoscalar scalar form fac- 
tor of the nucleon. Although not directly measurable this form factor <Tjv(t) is very interesting 
since it quantifies explicit chiral symmetry breaking effects in the nucleon (ajv(O) = mdM/dm 
is the celebrated pion-nucleon sigma-term). Moreover, its spectral function \ma N (t) is com- 
pletely dominated by the low-mass two-pion continuum due to the abovementioned threshold 
enhancement effect with no visible scalar resonance structures (see Fig. 4 in ref. [7]). 
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Fig.l: One-loop diagrams contributing to the imaginary parts of the isoscalar scalar and isovec- 
tor electromagnetic form factors of the nucleon. Dashed and solid lines denote pions and nu- 
cleons, respectively. The wiggly line represents the external scalar or vector source. The heavy 
dot symbolizes the second-order chiral nn NN -contact vertex proportional to the low-energy 
constants 01,2,3,4- The combinatoric factor of the last two diagrams is 1/2. 

Let us begin with recalling the basic definition of the isoscalar scalar form factor a N {t) of 
the nucleon: 

(N(p')\mqq\N(p)) = a N (t) u{p')u{p) . (1) 

Here, m = {m u + m^)/2 denotes the average light quark mass, t = (p' — p) 2 is the Lorentz- 
invariant squared momentum transfer and u{p) stands for a Dirac spinor. Similarly, the matrix 
element of the isovector vector current defines the isovector Dirac and Pauli form factors: 

(N(p')\qY^q\N(p))=u(p') 

where M = 938.27 MeV denotes the proton mass. In what follows, we will work with other 
linear combinations of F^ 2 (t), the isovector electric and magnetic (Sachs) form factors: 

G V E (t) = Fr(t) + ^ f2 F 2 v (t), G V E (0) = \, (3) 
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7 " fl(t)+ 2M (?, '"^ F2(t) 



T a u(p) , 
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GX,(t)=F?(t) + F?(t), 



M 



(0) 



= 2.353, 



(4) 



which obey at t — the normalization conditions written in eqs.(3,4). As usual, /i p = 2.793 
and n n = —1.913 denote the proton and neutron magnetic moments (in units of nuclear 
magnetons). In the complex t-plane the form factors cr N (t) and G^ M (t) have a right hand 
cut along the positive real axis starting at to — 4m^ related to the opening of the two-pion 
threshold. The imaginary parts Imo"Ar(t + i0 + ) and ImG^ M (t + i0 + ) of these nucleon form 
factors (equal to the half discontinuities across the cut) are also called their spectral functions. 

The contributions to the spectral functions lm a N(t) and lmG^ M (t) calculated up to two- 
loop order in chiral perturbation theory can be grouped into four different classes. The first 
class of contributions comes from the one-loop diagrams with leading order relativistic ttN- 
interaction vertices and propagators (represented by the first and second diagram in Fig. 1). 
For these one-loop graphs we can take the fully relativistic expressions given in ref.[5] and 
expand them in inverse powers of the large nucleon mass M up to order 0(M~ 2 ). This way 
one gets for the imaginary part of the isoscalar scalar form factor of the nucleon: 



Imcrjv(0 = 



3g 2 A m l 



(t - 2ml) A(t) 
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with the auxiliary functions: 
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t — Ami v 
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The same procedure gives for the imaginary parts of the isovector electric and magnetic form 
factors of the nucleon: 
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ImG^W = ^^ t {9AQ 2 m + ^[9l(^ml-At) + Q' 

+ ^ {3t2 ~ 12tml + 8m ' )A{t) }- (8) 

It is important to note that the expressions in eqs. (5,7,8) possess the correct threshold be- 
havior, namely lma N (t) ~ yjt — 4m 2 and lm.G\ M (t) ~ (t — 4m 2 ) 3 / 2 . In order to achieve 
this property in all three subleading term proportional to M 3 had to be kept in 

the isovector electric spectral function ImGg(t). One easily verifies that the coefficients of 



Q — yjt — 4m 2 appearing in the threshold expansion of ImG^ M (t) vanish to the respective 
order in the 1/M-expansion. The crucial point about the representations eqs. (5,7,8) is not to 
further expand the function A(t) in powers of 1/M as it is done (implicitly) in heavy baryon 
chiral perturbation theory [6]. In fact this function A(t) incorporates the anomalous (loga- 
rithmic) singularity at t c = 4m 2 — m^/M 2 , determined as that t- value where the argument of 
the arctangent-function becomes equal to ±i. As a matter of fact the 1/M-expansion of the 
function A(t) would destroy the correct analytical structure and moreover generate singular 
terms of the form Q~ n . From the point of view of chiral power counting the function A(t) 
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resumes an infinite string of terms starting at zeroth order. We note aside that the one-loop 
calculation of ref.[8] employing a so-called infrared regularization scheme gives naturally the 
correct threshold behavior of ImG^ M (t) since there the spectral functions are identical to 
those of the fully relativistic framework [5]. 

The next class of contributions is given by the one-loop diagrams with one second order 
chiral 7T7riViV-contact vertex proportional to the low-energy constants 01,2,3,4 [9]. We obtain 
from the last diagram in Fig. 1 the following contributions to the spectral functions: 



Im crjv(0 



ImGl 



mlQ 



Am\{c 2 + 3c 3 - 6ci) - t(c 2 + 6c 3 ) 
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c 4 Q 3 Vt 



lmG v 
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c 4 MQ 3 



(9) 



(10) 



192ttM/2' — mvv 4 87r j2^' 

where we have included also the first relativistic 1/M-correction (as far as it does not vanish). 
The separation of the 7r7riViV-contact vertex into an isoscalar part proportional to 0^2,3 an d 
an isovector part proportional to C4 is obvious from eqs.(9,10). The same C4-contribution to 
ImG^it) has been written down first in eq.(10) of ref.[6]. 
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Fig. 2: Two-loop diagrams contributing to the imaginary parts of the isoscalar scalar and isovec- 
tor electromagnetic form factors of the nucleon. The grey disk symbolizes all one-loop diagrams 
of elastic nN -scattering. The combinatoric factor of the first two diagrams is 1/2. The vertical 
dotted line indicates that cut which splits off the (one-loop) scalar or charge form factor of the 
pion from these factorizable diagrams. 

Now we turn to the two-loop diagrams. These are symbolically represented by the left graph 
in Fig. 2. The grey disk should be interpreted such that it includes all one-loop diagrams of 
elastic 7riV-scattering. Exploiting unitarity the imaginary parts \m.a^{t) and \mG\ M (t) can 
calculated from these two- loop diagrams as integrals of the product of scalar / vector source — > 
2n and 2n — > NN transition amplitude over the Lorentz-invariant two-pion phase space. 
When evaluated in the irir center-of-mass frame the pertinent two-pion phase space integrals 
become proportional to simple angular integrals j\ dx. Consider now that cut which splits 
a two-loop diagram into a tree-level source — > 2n coupling and a one-loop diagram of elastic 
7riV-scattering. Putting together all such terms we find the following two-loop contribution to 
the imaginary part of the scalar form factor of the nucleon: 



lma N (t) = 



3mz 



g\(l0ml -At) -ml 
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This expression is identical to the 27r-phase space integral 3m 2 (327^)^ J Q ff /2 dk Re g+(ik, t) 
with g + (u,t) the isospin-even non-spin-flip one-loop 7riV-amplitude [9, 10]. The analogous 
contribution to the isovector electric spectral function reads: 

lmG v (t) = AKiZnl-t) ^^^^ 



192tt 3 /M Jo " k 

1 i q\t , 9 , o^Q 4 , 9 . 16m^ A t 2 

x ln 4^ + <s^{ i< 257(2 + 26 ° 4 "* - 636S <> 

,9a, 
4 



10(Ji + J 2 )(2m 2 -i) - 3d 3 Q 2 + Wd 5 m 2 n . (12) 



l -Q 2 (13t - 48m 2 ) + ^(684m 4 + 73m 2 * - 16t 2 )l 
25 J 

+ <y 

480tt/|v^ 

The term in the last line originates from the third order scale-independent low-energy constants 
dj introduced in ref. [9] which subsume already chiral logarithms ln(mJA). The remainder in 
eq.(12) is identical to the 27r-phase space integral (Svrv^) -1 J®q/ 2 dkRe [— ikg~(ik, t)] with 
g~(u,t) the isospin-odd non-spin-flip one-loop 7riV-amplitude [9, 10]. The analogous contri- 
bution to the isovector magnetic spectral function reads: 

+ In + #(30m^ - 64mJ _ W (13) 

24^t ^ft-2m^ 24 v 71 1 J v y 

Again, this expression is identical to M(327rv^)^ 1 S-q/2 (Q 2 ~ 4/c 2 )Re h~(ik, t) with h~(uj, t) 
the isospin-odd spin-flip one-loop 7riV-amplitude [9, 10]. We emphasize the correct threshold 
behavior of the two-loop contributions: \vna N {t) ~ Q according to eq.(ll) and ImGg M (f) ~ 
Q 3 according to eqs. (12,13). In the latter two cases one easily verifies that the coefficients of 
Q vanish identically. 

The full set of two- loop graphs contains the two factorizable diagrams exhibited in Fig. 2. 
Both cuts of the intermediate pion-pair generate a contribution to the imaginary parts, accord- 
ing to the rule: Im^i^) = Re z\ • Im z 2 + Im z\ • Re z 2 . Only one term of this sum is accounted 
for in the expressions eqs. (11, 12, 13). The other term arises from that cut which splits off the 
one- loop scalar or charge form factor of the pion [11] (as symbolized by the vertical dotted 
line in Fig. 2). Naturally, the scalar form factor of the pion enters the scalar spectral function: 

(14) 

Here, (rg) n denotes the mean square scalar radius of the pion. The low-energy constant l 3 ~ 3 
[11] shows up in the one-loop representation of the pionic sigma-term: a n (0) = rhdm 2 /dm = 
m 2 + m\{\ — l 3 )/327r 2 f 2 ~ 0.99m 2 [11]. On the other hand side, the (one-loop) pion charge 
form factor [11] enters the isovector electric and magnetic spectral functions: 

lmGUt} = wm„yA q \- A{u ~ Sml) + Q2 1 " it" " 2mifA(t} } 

^N^'v™--^^}' (15) 
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ImG " (t) = w(4urA {27rf ' ] {rJ ' t + J ~ 2m ' ~ 27t "' ~^r\ • (16) 

where (r 2 h ) n denotes the mean square charge radius of the pion. In order to guarantee the 
correct threshold behavior ImGg(t) ~ Q 3 we have kept the first relativistic 1/M-correction 
(see eq.(7)) in the first factor of eq.(15). 
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Fig. 3: The spectral function Ima^{t) of the isoscalar scalar form factor of the nucleon multi- 
plied with ml/t 2 . The dashed- dotted line gives the one-loop result eq.(5). The dashed line 
includes in addition the C\^^-term eq.(9) and the full line includes furthermore the two- 
loop contributions eqs. (11,14). The upper dotted line shows the empirical spectral function 
milma N (t)/t 2 ofref.[7j. 

Let us add a remark on chiral power counting of the spectral functions. Counting the quan- 
tities y/i, m w and Q = \Jt — Am 2 as small momenta one deduces that two-loop contributions 
to Imajv(t), ImGg(f) and ImGj^t) are of fifth, fourth and third order in small momenta, 
respectively. In the latter two cases the electric/magnetic coupling of the (virtual) photon 
absorbs one/two powers of small momenta. 

Let us now turn to numerical results. We use the parameters: = 139.57 MeV (charged 
pion mass), f n = 92.4 MeV (pion decay constant) and qa = 1-3 (equivalent to a 7riV-coupling 
constant g n N = gAM/ f n = 13.2 [13]). For the second and third order low-energy constants q 
and dj we choose the average values of three fits to 7riV-phase shift solutions given in ref. [9]: 
c_i = -1.4 GeV -1 , c 2 = 3.2GeV-\ c 3 = -6.0 GeV" 1 , c 4 = 3.5 GeV" 1 , di + d 2 = 3.0 GeV" 2 , 
ds = —3.0 GeV -2 and d 5 = 0.1 GeV -2 . Most of these values are compatible with the low- 
energy constants q and dj obtained in ref. [12] in a one-loop analysis of 7riV-scattering inside 
the Mandelstam triangle. The mean square charge radius of the pion is accurately measured 
in elastic pion-electron scattering: (rl h ) 7T = (0.44 ± 0.01) fm 2 [14]. A value compatible with 
that has been found in ref. [15] in a two-loop analysis of the low-energy 7r~e~-scattering data. 



6 



An improved determination of the mean square scalar radius of the pion has recently been 
given in ref.[16]: (r|)?r = (0.61 ± 0.04) fm 2 . We use throughout the central values of these 
empirical pion mean square radii. 
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Fig. 4-' The spectral function ImG^(t) of the isovector electric form factor of the nucleon di- 
vided by t 2 . The dashed-dotted line gives the one-loop result eq.(7). The dashed line includes 
in addition the c^/M-term eq.(10) and the full line includes furthermore the two-loop contri- 
butions eqs. (12,15). The upper dotted line shows the empirical spectral function ImG^it) /t 2 
ofref.[2j. 

In Fig. 3 we show the spectral function Im<7jv(£) of the isoscalar scalar form factor of the 
nucleon multiplied with a weighting factor m\/t 2 . The dashed-dotted line gives the one-loop 
result eq.(5). The dashed line includes in addition the Ci, 2 ,3-term eq.(9) and the full line 
includes furthermore the two-loop contributions eqs. (11,14). The upper dotted line shows the 
empirical isoscalar scalar spectral function m\ lma^(t)/t 2 of ref. [7]. It has been determined 
on the basis of the 7T7t — > NN s-wave amplitude fl(t) [2] and 7T7r-scattering data tied together 
with Roy-equations etc. (see eq.(3) in ref. [7]). One observes a substantial improvement when 
including the next-to-leading order c 12 ,3-term and the two-loop contributions. The major 
effect comes evidently from the large isoscalar 7T7riViV-contact couplings, in particular from 
c 3 . The height of the peak at t ~ 5m 2 is however still underestimated by about 20% in two- 
loop approximation. Higher order 7T7r-rescattering effects etc. are necessary in order to close 
this remaining gap. Given the pattern in Fig. 3 one can expect significant effects (in the right 
direction) already from the two-loop diagrams with one vertex proportional to the (numerically 
large) second order low-energy constants 01,2,3,4. Note that a complete calculation of elastic 
7r./V-scattering to chiral order four which includes the pertinent one-loop diagrams with one 
Ci,2,3,4-vertex has been recently performed in ref. [17]. For comparison, similar deficiencies of 
the two-loop approximation of chiral perturbation theory have been observed in ref. [18] for 
the imaginary parts of the pion scalar and charge form factors. 
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Next, we show in Fig. 4 the spectral function ImGg(t) of the isovector electric form factor 
of the nucleon weighted with 1/t 2 . The dashed-dotted line gives the one-loop result eq.(7). 
The dashed line includes in addition the Czt/M-term in eq.(10) and the full line includes 
furthermore the two-loop contributions eqs.(12,15). The upper dotted line corresponds to 
the empirical isovector electric spectral function lm.G\(t)/t 2 of ref.[2]. Modulo a kinematical 
factor Q 3 /8My/i it is determined by the product of the nn — > NN p-wave amplitude f]_(t) [2] 
and the (time-like) pion charge form factor measured in the reaction e + e~ — > 7r + 7r~ (see eq.(7) 
in ref.[6]). In the case of the isovector electric spectral function the two-loop contribution 
is considerably larger than the c 4 /M-term (of the same chiral order). The peak at low itit- 
invariant masses, t ~ 7m 2 ,, gets effectively enhanced by a factor of about 1.2 by this additional 
ttN- and 7T7r-rescattering dynamics. The other prominent feature of the empirical isovector 
electric spectral function ImGg(t), namely its rise to the p(770)-resonance, is clearly absent 
in (chiral) perturbation theory. 
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Fig. 5: The spectral function ImG\j{t) of the isovector magnetic form factor of the nucleon 
divided by t 2 . The dashed-dotted line gives the one-loop result eq.(8). The dashed line includes 
in addition the c^-term eq.(10) and the full line includes furthermore the two-loop contributions 
eqs. (13,16). The upper dotted line shows the empirical spectral function lmG\it)jt 2 of ref.[2]. 

Furthermore, we show in Fig. 5 the spectral function ImG^(f) of the isovector magnetic 
form factor of the nucleon weighted with 1/t 2 . The dashed-dotted line gives the one-loop result 
eq.(8). The dashed line includes in addition the c^term in eq.(10) and the full line includes 
furthermore the two-loop contributions eqs. (13,16). The upper dotted line corresponds to the 
empirical isovector magnetic spectral function lmG^(t)/t 2 of ref.[2]. Modulo a kinematical 
factor Q 3 /8y / 2t it is determined by the product of the irir — > NN p-wave amplitude f*(t) 
and the (time-like) pion charge form factor (see again eq.(7) in ref. [6]). In the case of the 
isovector magnetic spectral function the two-loop contribution and the next-to-leading order 
C4-term are approximately equal in magnitude. Both together enhance the peak slightly above 
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threshold considerably by a factor of about 1.5. Again, the strong rise of the empirical isovector 
magnetic spectral function ImG)^(f) to the p(770)-resonance cannot be reproduced in (chiral) 
perturbation theory. 

For a more complete description of the dynamics governing the empirical isovector electro- 
magnetic spectral functions the introduction of an explicit p(770)-resonance is indispensable. 
For that reason we add a phenomenological p-meson exchange contribution of the simple form: 

Im G EM (t) - {m 2_ t)2 + tT 2 {t) , r,(t) - — , (17) 

which of course includes the proper energy-dependent p — > 2n decay width T p (t). We take 
the value m p = 769.3 MeV [19] for the p-meson mass and coupling constant g pn = 6.03 is 
determined from the empirical decay width T p (m 2 ) = 150.2 MeV [19]. Furthermore, b E = 
1.0 and hu = 3.5 are numerical parameters which we have adjusted to the height of the 
p(770)-resonance peak (at t ~ 28m 2 ) of the empirical isovector electric and magnetic spectral 
functions ImGg M (i)/f 2 . 
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Fig. 6: Spectral functions of the isovector electric and magnetic form factor of the nucleon 
weighted with 1/t 2 . Shown are ImG^(t)/t 2 (lower three lines) and ImG^it) /t 2 (upper three 
lines). The full lines give the two-loop results of chiral perturbation theory and the dashed lines 
include in addition the phenomenological p-meson contribution eq.(ll). The dotted curves 
correspond to the empirical spectral function ImG\ M {t)/t 2 of ref.[2j. 

Finally, the lower (upper) lines in Fig. 6 show the isovector electric (magnetic) spectral 
function weighted with 1/t 2 . The full lines give the respective two-loop result of chiral pertur- 
bation theory and the dashed lines include in addition the phenomenological p-meson exchange 
contribution eq.(17). One observes that the empirical isovector electric and magnetic spectral 
functions (dotted lines) can now be reproduced fairly well. The agreement is rather satisfac- 
tory for the electric one. In the case of the isovector magnetic spectral function \mG\jit) /t 2 
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some overshooting and undershooting of the empirical values below and above t ~ 13m 2 is 
visible in Fig. 6. This may point towards the relevance of higher order 7T7r-rescattering effects, 
or it may just result from the simplicity of the p(770)-resonance parameterization eq.(17). In 
any case one should not overinterpret the small differences between the dashed and dotted 
lines in Fig. 6 since there may be a potential danger of double-counting when adding together 
perturbative chiral pion loops and a finite width p(770)-meson exchange term (representing 
non-perturbative pion-dynamics). 

In summary, we have calculated in this work the imaginary parts of the isoscalar scalar and 
isovector electromagnetic form factors of the nucleon up to two-loop order in chiral pertur- 
bation theory. We have performed the non-relativistic 1/M-expansion such that the correct 
threshold behavior, lma N (t) ~ — 4m 2 anc ^ I m ^B,Af(0 ~ — 4m;;) 3 / 2 , is ensured. The 
two-loop contributions and the next-to-leading order Ci^^-terms magnify the strong enhance- 
ment slightly above the two-pion threshold. Higher order 7T7r-rescattering effects etc. are still 
necessary to close a remaining 20%-gap to the empirical isoscalar scalar spectral function 
Im<Tjv(t) [7]. After adding a phenomenological p(770)-meson exchange term to the respective 
two-loop results the empirical isovector electric and magnetic spectral functions ImG^ M (t) 
[2] can be reproduced fairly well. 

I thank Ulf-G. MeiBner for his tabulated values of the empirical isovector electromagnetic 
spectral functions \mG\ M {t) and for critical reading of the manuscript. 
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